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This article describes the Josephson tunneling from time-reversal symmetry-breaking states and 
compares it with that from time-reversal invariant states for both twinned and untwinned crystals 
and for both c-axis and basal-plane currents, in a model for orthorhombic YBCO. A macroscopic 
invariance group describing the superconducting state of a twinned crystal is introduced and shown 
to provide a useful framework for the discussion of the results for twinned crystals. In addition, a ring 
geometry, which allows s-wave and d x 2_ y 2-wa,ve superconductivity in a tetragonal superconductor 
to be distinguished on the basis of symmetry arguments only, is proposed and analyzed. Finally, 

, an appendix gives details of the experimental Josephson tunneling evidence for a superconducting 

■ state of orthorhombic ux 2 + vy 2 symmetry in YBCO. 
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I. INTRODUCTION 



Recently a number of Josephson tunneling cxpcriments0iH , i]i , il@ji'i , 0llllll@llilll have been performed on the high- 
Tc material YBCO to investigate the symmetry of its order parameter. Conventional superconductors all have order 
■ parameters with s-wave symmetry. On the other hand, many theoretical models of high temperature superconductivity 
based on a description of the pairing of the electrons in approximately tetragonal copper-oxide layers (such as occur 
in YBCO) find superconducting order parameters of symmetry d x 2_ y 2 (see Refs. 16| and |l7] for reviews). Thus, 
experiments have tended to focus on distinguishing between s and d x 2_ y 2 symmetries for the superconducting order 
parameter. 

It is of interest to note that this focus on the YBCO order parameter as a linear combination of d x 2_ y 2 and 
s-wave order parameters can also be arrived at from a study of existing experimental results only, without appeal 
to current microscopic theories. If the transition to the superconducting state is continuous, then according to the 
C| Landau theory of second-order phase transitions, the superconducting order parameter must transform according to 
one of the eight irreducible representations of the orthorhombic point group of YBCO. It is shown in Appendix A, 
that the observation of nonzero Josephson tunneling currents along all three principal crystallographic directions in 
YBCO (e.g. as in Refs. |l|, ||, |ll|) is sufficient to establish that the YBCO superconducting order parameter is a linear 
combination of d x 2_ y 2- and s-wave components. 

An important concept developed in this paper is that of the macroscopic symmetry (i.e. macroscopic invariance 
group) for the superconducting state of a twinned crystal. This macroscopic invariance group (which is different from 
the microscopic invariance group describing the superconducting symmetry of an untwinned single ctusta^Js. Ihc. 
symmetry group determined by Josephson tunneling experiments on twinned crystals. The majoritya'B'QBaQu'BEiJ 
of Josephson tunneling experiments on twinned YBCO crystals have determined the macroscopic symmetry of the 
superconducting state to be the same as that of a tetragonal d x 2_ y 2 superconductor. 

According to Ref. [II] (see also Ref. [ll]) , two types of twin boundaries are possible within a Ginzburg-Landau model 
of YBCO in which the individual twins are in d x 2_ y 2 + s states: the superconducting state can have either even 
or odd reflection symmetry with respect to the twin boundary. Ref. [l8| has shown that the experimental result of 
macroscopic d x 2_ y 2 superconductivity in twinned samples is equivalent to determining that the twin boundaries in 
YBCO have odd reflection symmetry. 

This leaves the question of the relationship between microscopic and macroscopic symmetry, which we discuss briefly 
below. As it turns out, a microscopic understanding of the energetics of the twin boundaries is necessary to determine 
whether or not macroscopic d x 2_ y 2 symmetry implies a microscopic symmetry which is predominantly d-wave (see 
also Ref. ^|). In this article, however, we prefer not to rely on microscopic arguments, and to consider all possibilities 
consistent with the requirements of symmetry. 

In principle, (e.g. assuming that difficulties concerning trapped flux can be satisfactorily eliminated) a direct ex- 
perimental determination of whether d x 2_ y 2 or s component of the order parameter is dominant can be made by 
performing a SQUID experiment on an untwinned single crystal of YBCO in the geometry suggested by Sigrist a: 
HiceE-3. Experiments on untwinned crystals have been performed by Wollman et al.uuQ and by Brawner and Ot 
and indicate that the microscopic symmetry of the order parameter is predominantly <i-wave. 
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The tricrystal ring experiments of Tsuei et alBB'a&}3& are of particular interest because of the precision with 
which problems associated with trapped flux can be eliminated, and because one of these experimentsEj has been 
carried out on Tl2Ba2CuOg+5, which is tetragonal and therefore not open to questions of interpretation related to 
the existence of twins. These experiments show convincingly rings containing a flux of one half an elementary flux 
quantum in their ground state; such rings must contain either one or three tt Josephson junctions. Based on simplified 
models for the angular dependence (on the angles between the Josephson junction and the crystallographic axes) of 
the Josephson current, a plausible argument is made that these results imply d x 2_ y 2 superconductivity. However, it 
would not contradict the basic principles of physics if, for the experimental geometries used up to the present, s-wave 
superconductivity were to give rise to a ring whose ground state contained half an elementary flux quantum. We 
therefore propose and analyze, in a section of our article below, a ring geometry where the observation of one half a 
flux quantum in the ground state can exclude, on the basis of experimental measurements and symmetry arguments 
only, the possibility that the superconducting state is s-wave. The case of d + is superconductivity (see below) is also 
discussed for this geometry. 

The experiments of Ref. [l^ do not appear to yield the d x 2 _ y 2 macroscopic symmetry obtained by other Josephson 
experiments performed on twinned crystals. A possible explanation for this apparent inconsistency has been given in 
Ref. ^ in terms of the existence of vortices in the weak links. Furt ier study of this experimental configuration has 
indeed found magnetic flux trapped in the grain boundary junctionsE3. 

A remaining puzzle among the Josephson results on YBCO-Pb junctions concerns the c-axis critical currents. 
Because of the orthorhombicity of YBCO, there should be a nonzero c-axis Josephson current from an untwinned single 
crystal even if the superconducting order parameter of YBCO is predominantly d x 2_ y 2-w&ve. In the measurements 



of Ref. 11, |12j, |13|, the c-axis current was found to be nonzero both for twinned and untwinned crystals of YBCO; 
the I c R n product (where I c is the critical current and S™ is the normal state resistance) is approximately an order of 
magnitude smaller for films than it isibr single crystalsE3, with the highest value being obtained for untwinned single 
crystalai3. On the other hand, effortsE-3 to detect Josephson tunneling in cases where the Josephson junctions have a 
relatively high value of R n have not been successful, perhaps due to thermal fluctuation effects. These varied results 
are not understood. According to Ref. H experimental measurements of basal-plane Josephson currents indicate 
that the twin boundaries have odd reflection symmetry; for such twin boundaries, the c-axis Josephson currents from 
twinned crystals will average to zero, assuming the two types of twins are present with equal weight. One possibility££l 
of accounting for a nonzero c-axis current is that the twins are not present with equal weight (perhaps as a result of 
strains); alternatives will be discussed below. 

It is perhaps useful to introduce the idea of "robustness" of the consequences of the macroscopic symmetry. The 
macroscopic invariance group of a twinned crystal is determined by assuming that all types of twins are present with 
equal weight. However, it may happen that, due to strain for example, imbalances occur where one type of twin has 
greater weight than another. A robust property is one which is not changed by a small imbalance of this type. It can 
be shown that the observation of a it phase shift in a basal-plane corner SQUID experiment on YBCO is a robust 
property, whereas the expected zero value of the c-axis Josephson current from a twinned sample having macroscopic 
d x 2_ y 2 symmetry is not. Thus, as noted in the previous paragraph, the c-axis results of Refs. [K5] could be 

explained in terms of an_unequal weighting of twin types. 

As noted by Li et al.c3 in a Ginzburg-Landau study of mixed s and d x 2_ y 2 superconductivity in an orthorhombic 
superconductor such as YBCO, the phase transition to the superconducting state will be, aspiming that this transition 
is continuous, to a time-reversal invariant state called the d+ s state. They also pointed outo that a second transition 
to a different superconducting state called the d + e^s state is possible at a lower temperature. This d + e l ^s state 
breaks time reversal symmetry. At present there is only one reporio of a second phase transition in YBCO. It is 
however possible, particularly in twinned samples, that such a second transition could be smeared out and not_easily 
visible. The d + e^s state is not stable in some microscopic theoriesEj; a recent weak-coupling calculation^, for 
example, finds that the d + s is always stable relative to the d + e^s state. Since there is at present no consensus 
concerning an appropriate microscopic theory of high temperature superconductivity, we do not exclude the d + e^s 
state from our considerations but explore its experimental consequences. 

This article begins by discussing the time- reversal symmetry-breaking d + e 1 ^ s state, the possible types of twin 
boundaries consistent with this state, and the macroscopic symmetries of twinned crystals having twins in the d+e^s 
state. Some results relevant to d + s states, which are a special case (<j) = 0) are included ^.completeness. _ . 

Broken time-reversal symmetry in YBCO has been searched for in optical experimentsoO, /xSR experiments^, 
and in Josephson tunneling experimentscl and has not been found. Nevertheless, because of the importance of the 
Josephson effect as a method of determining the order parameter symmetry, it is essential to have a clear understanding 
of the experimental manifestations of broken time-reversal symmetry for experiments performed on both untwinned. 
and twinned crystals. In, related work, some effects of broken time-reversal symmetry at crystal grain boundaries^ 
and at twin boundariescj have recently been considered by other authors. 

The Josephson current through a junction can be written as 
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/ = / c sin(0-0 c ), (1.1) 

where I c > is the Josephson critical current, and cj> is the gauge invariant difference between the phases of the super- 
conductors on the two sides of the junction. For normal junctions (e.g. junctions between conventions 
conductors) the characteristic phase of the junction can always be taken to be <p c — 0. ExperimentsB'E 
on high temperature superconductors on the other hand, give evidence for circuits containing an odd number of tt 
junctions, i.e. junctions having a characteristic phase of (p c — tt. The calculations performed below show, for both 
twinned and untwinned crystals, that a Josephson junction involving a superconductor in a state which breaks time- 
reversal symmetry can have a characteristic phase (pc which is neither zero nor tt. A similar conclusion concerning (p c 
was previously reached in Ref. Bl 



In sections [V] and VI , the Josephson currents for twinned samples are calculated by averaging over currents from the 
individual twins; the result of this averaging is shown to be consistent with the expectations based on the appropriate 
macroscopic invariance group for the twinned sample. 

In this article, it is assumed that the orthorhombic basal-plane unit cells of neighboring twins are at right angles 
to each other, cf. Fig. [|, whereas in fact they are rotated towards each other by a small angle. This misorientation 
was considered carefully in Ref. Il8] and can be shown not to alter the results for Josephson tunneling discussed here. 



II. TWIN BOUNDARIES BETWEEN STATES BREAKING TIME-REVERSAL SYMMETRY 

The Landau free energy density describing coupled s- and d x 2_ y 2-wave superconductivity in an orthorhombic 
superconductor such as YBCO is^J 

fh = a s \ip s \ 2 + a d \4> d \ 2 

+ (-l) e a(V^rf + 1>M + MPd + CVrf) 

+ {-i) e \jh\M 2 + foMiMd + 1>M, (2.1) 

where ip s and ipd denote the s- and d x 2_ y 2-w&ve components of the order parameter, respectively, where a s = 
a' s (T — Tgo), ad — a' d (T — Tdo), a' s ,a' d > 0, T is the temperature, and where e £ {1,2} indicates whic h of the two 
types of neighboring twins is described, cf. Fig. m The appearance of the (— l) e factors in Eq. (2.1), which is a 



generalization of the free energy presented in Ref. p3|, is related to the assumption that the twin boundary between 
two twins is a plane of reflection symmetry of the underlying crystal latticeta. With the definition of the coordinates 
in Fig. [I], the d x 2_ y 2 component of the order parameter, ipd, changes sign under reflection in the twinning plane. To 
ensure that the free energy of the reflected state is the same as that of the original state factors of (— l) e are included 
in Eq. (0). 

For a = fli = (3 2 = 0, with j3 a > 0, fld > 0, and (3 s fld — (Ai — 2|/?o|) 2 > for thermodynamic stability, the free energy 
( |2.l|) describes a tetragonal system undergoing a continuous phase transition at T = max^o, T s0 ) to either a d-wave 
(ipd 7^ 0, ip a = 0) or an s-wave (ip a ^ 0, ipd — 0) superconducting state. A second transition to a mixed (ipd ^ 0, 
ip s ^ 0) state occurs when a s f3d = «d(/?4 — 2\Po\) or adf3 s — a s (/?4 — 2|/3 |), respectively. We define the relative phase 
(p ad between ip s and ipd through 

ip d =\ipd\j* d , i>s = |V s |e^ +0 ' d) ■ (2.2) 

In the tetragonal case, <p s d is determined by the sign of {3q : 4>sd £ {0,7r} for (3q < 0, corresponding to a time-reversal 
invariant d + s state, (p sd € { 7r /2, 37r/2} for [3q > 0, corresponding to a time-reversal symmetry-breaking d + is state. 

The terms with coefficients a, (3\, and /?2 describe the coupling of the s- and d-wave components of the order pa- 
rameter which exists in an orthorhombic crystal. A continuous phase transition from the normal to a superconducting 
state occurs in the orthorhombic system when the smallest eigenvalue of the quadratic form in Eq. (^J) goes through 
zero. For a <C 1, i.e. weak orthorhombicity, and ad < a s the transition is at temperature Td ~ X^o + a 2 /a s a' d to 
a state whose d-wave component ipd grows with decreasing temperature approximately like \ipd\ 2 — ^'di^d — T)/ j3d, 
where [3d ~ [3d- The s-wave component ip s can be found as an expansion in powers of ipd, and to lowest order is 

ips = (~\y +1 (a/a s )ip d . (2.3) 

For this solution, ipd and ip s have the same phase (modulo tt) and the superconducting state is called a d + s state. 
For la/a^l <§; 1, \ip s \ \ipd\, and the superconducting state is said to be predominantly d-wave. For a s < ad, on the 
other hand, a transition occurs at T s w T s q + a 2 /ada' s to a predominantly s-wave superconducting state. 
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To determine the relative phase </> s d between ip s and 4>d, cf. Eq. (2.2), we set a > and, for simplicity, extend 
the ass ump tion of weak orthorhombicity by neglecting /3j and /?2 in comparison with /3q. By minimizing the free 
energy (2.1) we find, for /3 > 0, the following conditions on the relative phase 4> s d- 



sin(/) sd = 0, 



(2.4) 



cos0 sd = (-ir +1 (2.5) 
4/301^11^1 

Condition (|2.4|) yields a relative phase of (j> s d = for e = 1, 4> s d = tt for e = 2, and a time-reversal in vari ant d + s 



superconducting state. This state is realized near the transition from the normal state (see also Eq. (2.3)) and, as 
the temperature is low ered , remains stable until \ipd\ and \ip s \ have grown sufficiently to reduce the fraction on the 
right-hand side of Eq. ( [2.5|) to unity. At the temperature T = with 1^(^)11^(^)1 = a/4/?o a second transition 
to a d+e l ^s state takes place. For the e = 1 twin, for example, \4>sd\ increases continuously from zero for temperatures 
below Tj,. The values <p s d = i\4> s d\ correspond to two energetically equivalent states of twin e = 1 that are related to 
each other by the time-reversal operation T. The d + e^s state thus breaks time-reversal symmetry. 

In the case that /3q < 0, on the other hand, the minimum of the free energy ( |2.l| ) is always attained for cf> s d — for 
e = 1 and 4>sd — ft for e = 2, so that only one phase transition occurs and the superconducting state is always the 
time-reversal invariant state d + s. 

Our considerations so far have shown that the relative phase <fi s d between the components of the order parameter in 



a given twin is determined by the bulk free energy density (2.1). Contributions to the free energy due to the interaction 



between the order parameters on opposite sides of a twin boundary have to be considered to describe the relative 
phases of the order parameter components across a twin boundary. For a predominantly d-w&ve superconductor the 
interfaceJree energy per unit area associated with interactions between twin 1 and twin 2, cf. Fig. [j], is to lowest order 
given byx3 

ff 2 = B(r d ^d2 + ^dir d 2), (2.6) 

where ipdi and ipd2 denote the ci-wave components of the order parameter in twin 1 and 2, respectively. For time- 
reversal invariant d + s states, the free energy is minimized by tpdi — ipd2 {^di = — Vto) for B < (B > 0), which 
implies that the overall superconducting states have the form 

= (lpdl,lpsl,lpd2, *Ps2) = (4>d,4>s,1pd, -^s) , B<0, (2.7) 

and 



*g = (lpdl,1psX,1pd3,^s2) = {^d,i>s,-^dAs), B>0. (2.8) 

The states ^ g and \t„ satisfy a^ g = and <r^ u = — ^ u and are thus even and odd with respect to a reflection a 
in the twinning plane, cf. Fig. ^. An investigation of Josephson tunneling currents for time-reversal invariant states 
has shownEj that the experimental results for YBCO are consistent with odd symmetry under a, so that the d-w&ve 
component is the same on either side of the twin boundary. 

The assumption that ipdi = 4'd2 remains correct for the time- reversal symmetry-breaking d + e 1 ^ s states yields two 
possible states in twin 2 for a given state in twin 1: 

*i = e i **(hM,e , >- i |Va|), (2.9) 
*a = e i **(hM,-e < *-|Va|) ) (2-10) 
*2 = e^(|M-e-^ d |^l). (2-11) 

Since ^ = — cr^x, and since the state ^>2 is related to the time reversed state of twin 1, defined as = 
e l ^ d (\fpd\,e~' l ^'" 1 \tp s \), by ^ = — cr^i, it is clear that all states have the same bulk free energy density. Which 
of these two states of twin 2 prefers to be in contact with the given state of twin 1 is determined by the surface energy 
of interaction of twin 1 with twin 2. Note that if the pair (^>i,^2) is favored, the overall state is odd with respect 
to a reflection a in the twinning plane; this type of twin boundary will be called an odd-symmetry twin boundary. 
Otherwise, if the pair (^1,^2) is favored, the overall superconducting state is invariant, to within a phase factor, 
with respect to a reflection a followed by time reversal T; this second type of twin boundary will be called a Ta twin 
boundary. For a given material, only one type of twin boundary will be present; this will be the one which contributes 
the least to the free energy. 
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III. MACROSCOPIC SYMMETRY 



Since a twinned sample will have some regions in which the twinning planes are (110) planes, and other regions 
where the twinning planes are (110) planes, it is important to examine how the superconducting order is transferred 
from one region to another. Consider the case where the twinning planes are odd-symmetry reflection planes. Within 
a given twin, say twin 1, the sup erco nductivity is described by its order parameter ^i, which is a mixture of d x 2_ y 2- 
and s-wave components (cf. Eq. (2.9)). This order parameter is invariant with respect to the orthorhombic point group 
operations (rotations of 7r about the orthorhombic principal axes and reflections in planes normal to the principal 
axes) and hence, from the point of view of its point group symmetry, can be represented by a rectangle, as in Fig. [jj. 
(This same rectangle can also be used to indicate the relative orientation of the a and b axes in the different twins.) At 
an odd-symmetry twin boundary, not only does the rectangle representing the order parameter change its orientation 
on reflection across a twin boundary, but it changes sign. Thus, in Fig. ||, the alternating plus and minus signs in the 
rectangles represent the different signs of the order parameters in the different twins. Also, the interface between a 
region of (110) twin boundaries and a region of (110) twin boundaries, as determined by experimenttil, is shown in 
Fig. H It is of interest that this interface is itself largely a twin boundary, with small breaks in this boundary. Thus, 
the superconducting order parameter should be transmitted from a region of (110) twin boundaries to a region of 
(1 10) twin boundaries in exactly the same way as it is transmitted from one twin to another within a region of, say, 
(110) twin boundaries. In the case of Ta twin boundaries, the superconducting order is passed from twin to twin and 
between regions of differently oriented twin boundaries in a similar way. 



A. Macroscopic invariance group for d + e'^s states 

Figure ^ also gives a simple picture of the macroscopic symmetry of a highly twinned superconductor with odd- 
symmetry twin boundaries. Note that Fig. |^(b) is obtained from Fig. ||(a) by a rotation of about the c-axis. 
However, Fig. |3|(b) is also obtained from Fig. ||(a) by changing the signs of all of the order parameters, i.e. by 
multiplying the state of the entire system by minus one. Thus, a rotation of this twinned crystal by ^ir about 
the c-axis is equivalent to changing the sign of the overall order parameter. In other words, the superconducting 
state of the twinned crystal is invariant under the operation e l7r C*4. Similarly, the superconducting state of the 
twinned crystal can be shown to be invariant with respect to all of the operations of the direct product group 
D^\d 2 ) x /, where / is the group containing the identity and the inversion, and D^\d 2 ) is the combined group 
D { l ) {D 2 ) = (E,C 2 ,2U 2 ,2C 4 e t7T ,2U^e") in the notation of Ref. ||. The direct product group d£\d 2 ) x / will be 



said to describe the macroscopic symmetry of superconductivity in the twinned crystal. Note that this macroscopic 
symmetry group lacks time reversal; the time-reversal operation creates a different superconducting state with the 
same free energy. 

The macroscopic invariance group, D^(D 2 ) X J, of a twinned crystal with twins in d + e^s states and with odd 
reflection-symmetry twin boundaries is thus the same, except for the lack of the time-reversaLaperation, as the invari- 
ance group of the d x 2_ y 2 superconducting state of an (untwinncd) tetragonal single crystalEa. If the characteristics 
of the Josephson effect are determined by the macroscopic invariance group of the twinned crystal, then we expect 
the c-axis Josephson current to be zero, and an offset flux (cf. Eq. ( |6.2[ ) below) of i$o in a corner SQUID experi- 
ment. These characteristics of the Josephson effect for the type of twinned-crystal superconductivity described in this 
paragraph are confirmed by detailed twin averaging calculations in Sections [v| and VI . 



The discussion of the case where the superconducting twin boundaries are of type Ta can be carried out in a 



similar way. In this case, the basic states of the two types of twins are and ty 2 , as given by Eqs. ( p. 9] ) and (2.11). 
Detailed consideration of the effects of the elements of D4 shows that the macroscopic invariance group of this state 
of the twinned crystal is D { 4 3) (D 2 ) x /, where D { 4 3) (D 2 ) = (E,C 2 ,2U 2 ,2Te ln C 4: ,2Te l7V U 2 '). Thus, the macroscopic 
symmetry here is the same as that of the d + is superconducting state of an (untwinncd) tetragonal crystal. Again, the 
time-reversal operation T by itself does not leave the state invariant, but rather creates a different state with the same 
free energy. The macroscopic symmetry here allows a nonzero c-axis Josephson current and an offset flux differing 
from i$o m a corner SQUID experiment; detailed calculations performed below indeed confirm these properties for 
a twinned crystal with Ta twin boundaries. 



B. Macroscopic invariance group for d + s states 

A previous papeiEl which presented detailed calculations of the Josephson currents from an orthorhombically 
twinned crystal, with twins in the time-reversal invariant d + s state, to a conventional superconductor, did not 
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exploit the idea of the macroscopic symmetry of the twinned cryst al. Two types of twin boun daries are possible, 



giving rise to superconducting states that are either odd, cf. Eq. ( |2.7j ), or even, cf. Eq. with respect to a 



reflection in the twinning plane. The case of odd reflection-symmetry twin boundaries is the same as the case of odd 
reflection-symmetry twin boundaries studied above, except that 4>sd — since the twins are in time-reversal invariant 
states. Hence, the time-reversal operation T is added to the invariance group, which becomes D^p (D2) x / x i?, where 
R contains the identity and time-reversal; this is also the invariance group for d x i_ y i superconductivity in a tetragonal 
superconductor. For the case of even reflection-symmetry twin boundaries, which can be represented schematically 
by assigning plus signs to all of the rectangles in Fig. |^, the macroscopic invariance group is D4 x I x R, which is 
also the invariance group for s-wave superconductivity in a tetragonal superconductor. Thus, in orthorhombically 
twinned crystals, with the twins in d + s superconducting states, odd reflection symmetry twin boundaries give a 
macroscopic symmetry identical to that of a tetragonal rf 2 .2_ y 2-wave superconductor, whereas even reflection-symmetry 
twin boundaries give a macroscopic symmetry identical to that of a tetragonal s-wave superconductor. Clearly, as 
emphasized in Ref. Q the qualitative characteristics of the Josephson effect in twinned crystals, which are determined 
by the macroscopic symmetry of the superconducting state, are a reflection of the symmetry of the twin boundaries, 
and have nothing to do with whether the basic d + s state of the twins is predominantly d x 2_ y 2--w&ve (\tpd\ ^ IV'sl) or 
predominantly s-wave (\ips \ 2> iV'dl)- 



C. Relation between twin-boundary symmetry and microscopic symmetry 

Given that Josephson experiments on twinned crystals of YBCO determine directly the twin boundary symmetry, 
it is of interest to ask to what extent this result can be used to infer the dominant microscopic order-parameter 
symmetry. 

As note d ab ove, if a phenomenological Ginzburg-Landau type of model is adopted for the twin boundaries (e.g. 
as in Eq. ( |2.6| ) above or as in Ref. |l9| ) they can have either even or odd reflection symmetry. Which case occurs 
depends on the signs and magnitudes of certain parameters in the twin-boundary free energy. The fact that the 
microscopic order parameter may be either predominantly d x 2_ y 2-w&ve or predominantly s-wave places no restrictions 
on these parameters at the phenomenological level. Thus, within the framework of Ginzburg-Landau type models, 
the identification of the twin-boundary symmetry does not lead to information about the dominant order parameter 
symmetry. 

Now consider the odd-symmetry and even-symmetry twin boundaries as shown in Fig. |2|. For the odd-symmetry twin 
boundary, the s-wave component of the order-parameter changes sign, whereas the d-wave component is continuous. 
On the other hand, for an even-symmetry twin boundary, the d-wave component of the order parameter changes sign, 
while the s-wave component is continuous. If it is assumed that either the d x 2_ y 2- or s-wave component of the order 
parameter is completely dominant, and that the lowest free energy state is achieved by having an order parameter 
that is continuous and varies as little as possible at the twin boundary, then predominantly d x 2_ y 2-w&ve (s-wave) 
superconductivity would lead to odd-symmetry (even symmetry) twin boundaries. The argument that the dominant 
order parameter should be as continuous as possible seems reasonable, but there are at least hypothetical situations 
(see e.g. Ref. ^) for which-it is not correct. .— . 

In YBCO, both infraredEi and Josephson tunnelingt^J measurements show a relatively strong anisotropy in the ab- 
plane penetration depth. Furthermore, a recent microscopic analysisEj of the penetration depth data indicates that the 
CuO chains play a significant role in the superconductivity. This suggests that neither the d x 2_ y 2 nor the s component 
of the order parameter is completely dominant, but that both have significant amplitudes. For an odd-symmetry twin 
boundary, the assumption that the order parameters vary continuously across the twin boundary implies that the 
c? 2 .2„j / 2-wave component of the order parameter remains approximately uniform across the twin boundary. The s-wave 
component, however, will have its bulk value ip s at distances greater than — £ s from the twin boundary (£ s is the 
s-wave coherence length), will go continuously to zero and change sign at the twin boundary, and will tend to the 
value — ips at distances greater than +ip s from the twin boundaries. 

If the condensation energy per unit volume gained by having a nonzero s-wave order parameter is e s , then the free 
energy per unit area of an odd reflection-symmetry twin boundary will, in this model, be approximately 

/odd = 2e s £ s + Sfodd- (3-1) 

Here a quantity <5/ dd is assumed to include all corrections necessary to arrive at the correct answer; these corrections 
may be due to interactions between the s and d-wave components of the order parameter, effects of crystal structure 
and doping level at the twin boundary giving condensation energies not equal to the bulk value, the effect of disorder 
at the twin boundary which might reduce the amplitude of the <i-wave component of the order parameter, etc. 
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Similarly, in the case of an even reflection-symmetry twin boundary, the twin-boundary free energy per unit area for 
our simplified model can be written 

/even = 2<€d^d ~t~ ^/cvcn- (^-2) 

Within this basic model of the energetics of the twin boundary, a knowledge of the relative condensation energies, 
and of the relative coherence lengths of the s- and d x 2_ y 2-w&ve components of the order parameter is required to 
determine whether the odd-symmetry or even-symmetry twin boundary is stable for a given relative weight of s and 
d x 2_ y 2 to the order parameter. We are not aware of the existence of good quantitative estimates of these quantities. 

In summary, our view is that, while it may be possible to make appropriate assumptions which would suggest 
that the observation of odd-symmetry twin boundaries implies a dominant d x 2_ y 2 component of the order parameter, 
the arffijmjwait is not conclusive. There is thus an important, role for experiments on untwinned single crystals of 
YBCOH'cHj'EI, or on tetragonal crystals such as Tl2Ba2Cu06+5c2l, since in these cases the complications due to twinning 
are eliminated. 

Finally, it should be noted that related issues have also been discussed in Ref. ||, where a somewhat different point 
of view is taken. 



IV. JOSEPHSON CURRENTS 



This section establishes a relationship between the-Josephson current across a junction and the phase derivative of 
the free energy by extending the approach of AndreeuLj (see Ref. [I?] for a review) to the case of a multicomponent order 
parameter. In order to calculate currents in the superconductor, gradient terms compatible with the orthorhombic 
symmetry of the system are added to the free energy density To lowest order, these are given by 



a X.jj, 



(4.1) 



where a € {x,y,z} and A, \x £ {d, s}. The coefficients of the gradient terms are real, satisfy 
components p a of the "momentum" operator are given by 



m° A , and the 



Pa 



(4.2) 



where A a is a component of the vector potential, c is the speed of light, and |e*| is the absolute value of the charge 
of a Cooper pair. The components j a of the current are obtained from the variation of F gla d with respect to A a , and 
are given by 



ja = X 2m 5 ~^ a ^" + C ' C ' ' 



(4.3) 



To derive a Josephson current-phase relation for our unconventional superconductor, we ge ner alize definition (|2.2|) 



by writing ip^ = \ifj tl \e l ^ d+<t '>* d \ where (pdd = 0, and 4> s d is the relative phase defined in Eq. (2.2). The phase 4>d can 
vary with position but, as pointed out in the previous section, 
assumption that \ip^\ is independent of position, we find: 



b s d is determined by the bulk free energy. Under the 



d(pd 2tt 
dx a $ 



A 



he 



(4.4) 



where 



file* 



rn 



A, fj. 



■|^A||Vvl C0S (^V"i ~ fad) ■ 



(4.5) 



The Josephson current across an interface between YBCO and a conventional superconductor like lead depends 
on the interface free energy associated with interactions between the s- and d-wave components ipuiP £ { s j^} °f 
the order parameter on the YBCO side and the conventional (s-wave) order parameter ippy, on the lead side. This 
interface free energy per unit area can be written as: 
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// = E G M(n)(^Pb+ c.c. 



(4.6) 



where G M (n) is real and, in general, depends on the direction of the surface normal n, relative to the YBCO crystal- 
lographic axes. There are also contributions to the interface free energy that depend on ip^ and 0pb separately. Since 
these terms do not contribute to the tunneling we neglect them here. 

The minimization of the total interface free energy with respect to variations in tp^ yields 



E" tt 2 I a CP«Vv) = -GA(n)^ Pb . 



(4.7) 



Insertion of Eq. fl4.7[ ) into (4.3) gives, for the Josephson current j across an interface between YBCO and lead with 
normal vector n, 



n J = 



A 



Pb 



(4.8) 



As in the case of conventional superconductivity, the current j is related to the derivative of the interface free energy 
fx with respect to the phase difference between the order pa ramet ers o n both sides of the junction. Substituting 
V-"Pb = \''pPb\e' l ^ Fb , and ipd and ip a from Eq. (2.2) into Eqs. (4.6) and (4.8), and setting <fi — <pd — <ppb yields 



36 2irc 



n J 



(4.9) 



V. c-AXIS TUNNELING 



Consider as a first example Josephson tunneling along the c-axis between YBCO and lead. The free energy per unit 
area // describing the interaction between the superconducting order parameters at the interface, which is normal to 
the c-axis, can be written as: 



(5.1) 

|). The coefficient 



where the constants (— l) e+1 D and S are appropriate realizations of the coefficients G M (n) in Eq. 
D is nonzero because of the orthorhombicity of YBCO. 

We start with the interface energy fz t i for twin 1 for the case that YBCO is in the time-reversal symmetry-breaking 
d + e^s state. Let twin 1 be in the superconducting state \&i defined in Eq. (2.9) and let ipp^ = |i/>pb|e l( ^ pb . With the 
definitions 



d= 2-D|^d||V>pb| . s = 2S|V> s ||V>Pb| , 6 = 6^-6 



Pb : 



we find 



with 



-r% cos ( 



l). 



(5.2) 



(5.3) 



sin 0i 



n = \/(d + s cos 6 ad ) 2 + (s sin 6 sd ) 2 
-ssm6sd/ri , cos 6% = (d+ scos6 s d)/ri 



(5.4) 



Note that for the time-reversal symmetry-breaking state (i.e. 6 s d ^ 0), the characteristic phase of the junction, 6\, 
will be neither nor n. 

If twin 2 is in the state ^2 = —0^ 1, c f. Eq. ( 2.10 ), so that the twin boundary is of the odd symmetry type, then 
the free energy //^ follows from Eq. fl5.3j) for fj t \ by replacing 6sd — » 4>sd + tt and d — > — d with the result: 



1 1.2 = -n cos ( 



- tt) = -/, 



1,1 



(5.5) 



Henc e, t he average of the interface energy over the two types of twins (fi,i + fi,2)/2 vanishes which implies, cf. 
Eq. (4.9), that the Josephson currents also average to zero. This result of the Josephson current averaging to zero for 
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odd symmetry twin boundaries is independent of the value of (f) s d and is thus the same for the time-reversal invariant 
and the time-reversal symmetry-breaking states. Hence, the result obtained above from the ideas of the macroscopic 
symmetry of a twinned crystal is confirmed: the element e I7r C4 of the macroscopic invariance group will not allow 
c-axis Josephson currents. 



Now suppose that twin 2 is in the state ^ = ~ <T _ ^ii cf- Eq. (2.11), so that the twin boundaries are of the Ta type 



The free energy jj i is then obtained from Eq. ( p.3|) with the substitutions <p s d —* —<t>sd + f and d — > —d, which yields 

/r,2 = -n cos (<j> + 4>i - 7r) . (5.6) 
The average of the interface free energies over the twins is then given by 

// = ^(fi.i + fi.2) = -risin^isin^. (5.7) 

Together with Eq. flt.SQ this yields for the current, 

n j = ^^ Cri sin 0i cos (j) ■ (5-8) 

Hence, there is a partial cancellation of the Josephson currents upon averaging over twins such that the critical 
current j c is reduced from the value for the untwinned crystal by a factor | sin </>i | . This result of partial but not 
total cancellation of the Josephson currents in twinned crystals resembles to some extent what has been observed 
in experiments, namely, that the highest values of the I c R n product have been observ ed , in untwinned crystalsEij, 
and that smaller values (which are however nonzero) are observed in twinned crystalstJ'EHrlla (see however Ref. fl4| ). 
Note finally, that the macroscopic invariance group found above for the case of Ta twin boundaries allows a nonzero 
Josephson current, in agreement with what has been found here. 



VI. BASAL PLANE JOSEPHSON CURRENTS 



tr 



erimental study of the basal plane Josephson currents is often carried out using the corner SQUID geome- 
ilTof Fig.g The total current /total, which is the sum of the currents through junctions A and B is 



/total = Ia sin(<^4 - 4> c .a) + Ib sin(0 B - </> c b) 



(6.1) 



The difference in the gauge invariant phase differences for the two junctions is <\>a — <t>B = 27r<E>/ <&0j where $ is the 
flux through the SQUID and $o is the elementary flux quantum. The average phase (4>a + <Ab)/2 is a free variable 
which adjusts itself so that the total current /total equals the current fed into the SQUID. The maximum value of 
/total with respect to variations of {4>a + <Ab)/2 is 



{I A + Ib) 



(1 



e 2 W 



Tj-(*-$of&et) 



1/2 



(6.2) 



where 



<JcA 



»cB 



offset 



2- 



(6.3) 



and 



e = {I A -I B )/{I A +I B ). (6.4) 

Note that / ma x is periodic in the flux $ with period <f>o- The quantity $ ffsot is the flux at which / max has its largest 
value. For tjwftjiormal junctions $ fi'sct = 0, whereas for a SQUID with one normal and one 7r junction, $ ffsot has the 
value $o/2EJaO. It is now shown that ^offset can have other values when the superconducting state is not time- reversal 
invariant. 

Consider a plane Josephson junction between twin 1 of YBCO and an s-wave superconductor such as lead. The 
junction is normal to the YBCO basal plane and its free energy fi\ per unit area is 

f IA (9) = - [D(9W dl + S(9)p Bl ] yj Pb + c.c. , (6.5) 
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where 8 is the angle between the normal to the interface n and the basis vector ai. By reflecting the state (ipdi, "4>si) 
of twin 1 with interface angle 8 in the twin boundary, a state (ipd2 , ?Ps2 ) = ( — ipdit^si) for twin 2 is obtained which 
has interface angle — 9 + tt/2 and the same free energy. Hence the interface free energy for twin 2 is given by 



Orthorhombicity implies that 



/i >( 2)(0) = - \-d(-9 + *-w d2 + s(~9 + 

+ c.c. . 



D(8 + tt) = D(-8) = D(8) , 
S(8 + tt) = S{-9) = S{9) . 



Pb 



(6.6) 



(6.7) 



With the aid of definition (2.9) for the superconducting state, the interface free energy for twin 1 can be written as 

//,i(0) = -r(0) 008(0-^(0)), (6.8) 

where cf> = cf>d — </>pb, and where 



r{9) 



d 2 (8) + P(8) 



m = 



r(0)' 



(6.9) 



with 



d{9) = 2(D(8)\^ d | +S(6)\tp s \ cos fa) |V>p b I 
s(8) = 2S(8)\i> s \\i>Vb\s^fa. 



(6.10) 



Com binin g Eqs. ([1.9]) and (6 



gives, for an untwinned single crystal, a Jos eph son current in the standard form 
of Eq. (1.1), with the characteristic phase of the junction 4> c (8) given by Eq. (6.9). For the time-reversal invariant 
d + s state, <j) s d = and the characteristic phase of the junction will be either zero or tt. Results appropriate to 
a tetragonal superconductor in a d x 2_ y 2--wave state are obtained from the above formula by taking S(9) — and 
D(8 + ^7r) = —D(8). This gives <j> c (8 + hw) = 4> c (9) + tt and an offset flux in a corner SQUID experiment of 
^offset = 5^0- F° r an orthorhombic crystal such as YBCO, there is however no symmetry- required relationship 
between D(9+ ^tt) and D(9), and also S(9) is nonzero in general. Nevertheless, provided the quantities d(9+ ^tt) and 
d(9) have appropriate signs, the d + s state (cf) s d = 0) will yield the value <I> offsot = m a corner SQUID experiment. 
Under these conditions, the superconducting state of orthorhombic YBCO may be said to exhibit predominantly 
d x 2_ y 2-wa,ve behavior. For the time-reversal symmetry-breaking d + e 1 ^ s state, <fi s d ^ 0, and 4> c (8) and 4> c (8 + ^tt) 
differ from zero and tt, and ^offset differs from zero or If \' l Ps\ <C \ip d \ and/or (f> s d <C ^tt, this difference may be 

small and difficult to detect experimentally. 

To describe basal-plane tunneling between a twinned sample of YBCO and lead, the interface free energy // is 
averaged over the two ty pes o f twins. The interface free energy fr?.( 8) of twin 2, cf. Eq. (3.6), with the superconducting 
state <J/2j given by Eq. ( 2.11 ), can be expressed in the form of Eq. (|6.8| ): 



/,2 



-r(9 + tt/2) cos(tt + cf> + <p c (9 + tt/2)) , 



which yields for the average fi(9) = (fi,x + fi.2~)/2 

/j(0) = -T*(0)cosfa-&(0)) ; 

where 



(6.11) 



(6.12) 



r- 2 (9) -- 

COS (j>2(9) 



d 2 d (0)+PM 
d d (9) 



r- 2 (6) 



, sin $2(8) 



Ss(9) 



(6.13) 



Here the ^-dependent functions have been separated into even (s-wave like) and odd (ci-wave like) terms under 
9 -► 8 + tt/2 by 
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x{6) = ~ [x{9) + x{6 + tt/2)] + 1 [x(6) - x{6 + tt/2)] 

= x s (0) + z d (0), (6.14) 

where a;(0) stands for or s(0). Also, since the two twins have been taken to be in states Vfi and W?, th e twin 



boundaries are what have been called above Ter twin boundaries. The above results (in particular Eq. (3.13)) show 
that for a twinned sample with Ter twin boundaries the characteristic phases of the Josephson junctions satisfy 
4>2{0 + \t) = —4>2{Q) + it- This yields an offset flux in a corner SQUID experiment of 

^offset = — *o + • (6.15) 

TT Z 

The macroscopic invariance group of a crystal with Ter twin boundaries, which is the same as that of a tetragonal 



crystal in a d + is state, does not constrain the offset f ux to be or ^tt, in agreement with Eq. (3.15). A value of 
^"offset = (0.98 ± 0.05)<f'o/2 is obtained in experimentsH on YBCO. Thus, either the YBCO twins are not in states 
of broken time- reversal symmetry (i.e. <^ s( j = 0), or the time-reversal symmetry-breaking is small (more explicitly, 
s s <C dd), or the twin boundaries are not of the type Ter, or time-reversal domains (see below) play a role. 

Alternatively, if the twin boundaries are of the odd reflection symmetry type, then twin 2 is in state given by 



Eq. ( |2.10D , and Eq. fl6.13| ) holds but with s s (8) replaced by s d (6). This yields <j> 2 {6 + |tt) = <j> 2 (9) + tt for the critical 
phase of the Josephson junction, and $ ffset = h^O- Thus, for odd-symmetry twin boundaries, the same value of the 
offset flux is obtained for the time-reversal invariant d + s state and for the time-reversal symmetry-breaking d + e'^s 
state. It is therefore not possible to determine whether or not YBCO is in a time-reversal symmetry-breaking state 
from the experimental results of Ref. ^ alone. These results can again be qualitatively understood in terms of the 
macroscopic invariance group of the twinned crystal, which for this case is the same, except for the absence of time 
reversal, as that of the d x 2_ y 2 superconducting state of a twinned crystal. 

There remains the question of the existence of domains which are time-reversal images of one another. Given 
a time-reversal symmetry-breaking state of a twinned crystal, a different state of the crystal having the same free 
energy can be obtained by applying the time-reversal operation. The various equations describing this state can be 
obtained from those presented in this and the preceding Section by reversing the sign of 4>sd- It may be that these two 
different states exist in different regions of the same crystal, in which case the crystal is said to contain time-reversal 
domains. If these time-reversal domains are sufficiently small, so that it is necessary to average over them to obtain 
the Josephson currents through a junction, then the behavior characteristic of broken time-reversal symmetry will 
disappear. This makes sense since a crystal containing a random distribution of domains which are time-reversed 
images of each other is macroscopically invariant under time reversal. Thus, the c-axis Josephson current found in 



the case of Ter twin boundaries, cf. Eq. (5.8), vanishes when averaged over time-reversal domains (since (f>i changes 
sign in going from a given domain to its time reversed image). Similarly, an average of the basal plane Josephson 
currents over time-reversal domains will yield junctions whose characteristic phases can only be either or tt. 



VII. PROPOSED GEOMETRY FOR RING EXPERIMENT 



Recentl^^njimber of order parameter symmetry determinations have been carried out using tricrystal ring mag- 
netometryuBaHEl These experiments have measured the flux through a superconducting ring containing three 
Josephson junctions. The observation of a spontaneous magnetization of the ring corresponding to a magnetic flux 
through the ring of ^$01 where $0 is the elementary flux quantum, shows that one or perhaps three of the Josephson 
junctions have a characteristic phase of tt. These experiments are impressive in eliminating the possible influence of 
spurious trapped flux, and in the precision with which the spontaneous flux is measured. 

The identification of a Josephson junction having a characteristic phase of tt was taken to be evidence for d x 2_ y 2 
symmetry. This conclusion was based on a simplified mode£3 of the angular dependence (on the angles between the 
Josephson junction and the crystallographic axes) of the Josephson current expected for d x 2_ y 2 superconductivity. 
The use of this simplified model is somewhat analogous to the description of an electronic band structure in a tight 
binding model restricted to nearest neighbor interactions only. Such a model is a good first guess, but its reliability 
is in fact unknown. A complete parametrization of the Josephson current is given in Appendix B, where it is seen 
that, even to "lowest order", the parametrization of Ref. |37] does not retain all terms. It should also be stated that, 
in its study of YBCO, Ref. || examined three different geometries all of which provided results consistent with the 
simplified d-wave model, and inconsistent with other possibilities which they raised. For the geometries studied in 
Refs. ^, [| [n], |2^, however, it would be extremely difficult to reach definite conclusions concerning the symmetry if 
a general model (similar to that developed in Appendix B) for the Josephson current were adopted. For this reason, 
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we propose below a different ring geometry which can distinguish between d x 2_ y 2 and s wave superconductivity on 
the basis of symmetry arguments only, and without reference to simplified models of the Josephson junction. 

It should also be noted that, as was the case with the corner SQUID experiments discussed above, the symmetry 
which is determined in these ring experiments on twinned crystals of YBCO is the macroscopic symmetry. Thus, a 
determination that the macroscopic symmetry is d x 2_ y 2 would lead directly to the conclusion that the twin boundaries 
have odd reflection symmetry, while further arguments are needed to determine the microscopic symmetry of the order 



parameter, cf. Section [II C. Experiments on tetragonal crystals, such as those of Ref. |20| on Tl2Ba2CuC>6+,5 are thus 
of special importance in that the complicating effect of twinning is eliminated. 

In view of the above, it would be useful to have an experimental test for s-wave versus d x 2_ y 2-wave superconductivity 
in tetragonal crystals which is dependent on symmetry arguments only. The experimental geometry of Fig. provides 
such a test. The superconductors labeled 5*1 and S2 are the same (e.g. Tl 2 Ba 2 Cu0 6+ 5), but 5*2 has its a and b axes 
oriented along the x and y axes of the figure, whereas the a and b axes of SI are rotated from x and y through an 
angle of 7r/4 about the c axis. 

Considerations of gauge invariance and time-reversal symmetry show that the sum of the free energies of the 
Josephson junctions A and B (to terms quadratic in the order parameter) has the form 

F = C A {^IA$>*2A + *1A*2A) + C B (^1B% B + *1B*2b) , (7.1) 

where "3/ia is the order parameter of superconductor 5*1 at Josephson junction A, etc. , and where ca and cb are real. 
If the superconductivity is s-wave, the invariance of the free energy with respect to rotations of the ring by n about 
the x axis yields the result ca = Cb; in this case both junctions are either it junctions or normal junctions, and the 
ring will not exhibit a spontaneous magnetic moment in its ground state. On the other hand, if the superconductivity 
is d-wave, the same symmetry argument requires that ca = —cb) hence in this case one junction is normal while the 
other is a tt junction and the ring will exhibit a spontaneous moment in its ground state. The flux associated with 
this moment is i$ i n the limit (LI c /c) ^> &o/(2n), where I c is the critical current and L is the self-inductance of 
the ring0S§. 

Note that in the geometry of Fig. || both a — and a — ir/4 yield ca = cb = for d x 2_ y 2-w&ve superconductivity 
so that these angles would not be appropriate choices for the experiment. Hence, values of a such as a ~ 7r/8 should 
allow a clear distinction to be made, on symmetry grounds only, between s-wave and d x 2_ y 2-w&ve superconductivity. 

In addition, it has been pointed out to usE3 that the model of Ref. || for the strong disordered limit of the Josephson 
junction yields ca = cb = for all a; thus, it would be prudent to have the least possible disorder in the junction. 

We now give a detailed analysis which includes three possible types of superconductivity for a tetragonal system, 
cf. Section O, namely d x a_ y a, s and d + is. Assume that in general the superconductivity can be described in terms of 
two complex order parameters, ipd and ip s , of d x 2_ y 2 and s symmetry, respectively. Then the free energy of junction 
A has the form 

Fa = G d 1p dlA 1pd2A + G s ^ slA 1p* s 2A + G sd ll) s iA^*d2A + Gds^dlA^tiA + cx - ' ( 7 - 2 ) 

where ipdiA is the order parameter of d x 2_ y 2 symmetry in superconductor SI at junction A, etc. , and where all 
coefficients G are real. The free energy for junction B is the same, except that Gd — > — Gd, Gds — ► ~Gd s , and the 
subscript A is replaced by subscript B everywhere. For the d + is state, the order parameter at junction A takes the 
form 

^diA = |V>dlA|e # ^ ,^ d2A = \ipd2A\e l * d2A , 
Aia = \AiA\e^ dlA+ ^ , yj s2A = A|^2A|e iW ^+^ ) , (7.3) 

where A = ±1 accounts for the freedom to choose either of the two d ± is states in superconductor S2. The order 
parameter at junction B is similarly expressed. Further analysis shows that the sum of the free energies of junctions 
A and B can be written as 

F = ~y~c Ia cos ^ a " ^cA) - ~^~Jb cos{(f>B - c b) , (7.4) 

where I a, Ib are positive. For pure d x 2_ y 2--w&ve superconductivity, 4> c a — 4> c b = ±t- For pure s-wave superconduc- 
tivity, tficA — 4> cB = 0. For d + is superconductivity, 4> c a and 4> c b are (as a result of nonzero G s d and Gds) different 
from zero (mod 7r). The quantities 4>a = 4>diA — <Pd2A and 4>b = 4>diB — 4 l d2B, when generalized to the usual gauge 
invariant form, yield 



-(^-«fe) = — 



(7.5) 
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where $ is the flux through the ring. The quantity \{cj>A + </>b) is at present undetermined; minimizing F with respect 
to it yields 



F = 



$0 (T 



e 2 + (l-e 2 )cos 2 



1 

2? 



where 



-j-($-$off S et) 



1/2 



(7.6) 



offset 



VcA 



2tt 



(7.7) 



and 



\Ia-Ib\ 
Ia+Ib 



(7.; 



Here, a quantity L/ 2 / (2c 2 ), representing the magnetic field energy associated with the current I circulating in the ring, 
has been added to F. The flux through the ring due to this current is = LI/c, where L is the ring self-inductance, 
and the total flux through the ring is $ = + <I> CX , where $ e x is the flux due to external sources. 

The ground state free energy of the ring in zero external flux, $ cx = 0, is found by minimizing F with respect to 
For sufficiently large values of the parameter L(Ja + /s)27r/(c c E>o), the term in LI 2 can be treated in perturbation 
theory. Also, without loss of generality, we can assume < (4> c a ~ 4>cb) < 7r. From these assumptions it follows that 
the flux through the ring in its ground state is $ ffsct- Thus, the flux through the ring in its ground state is zero 
for an s-wave superconductor, for a d x 2_ v 2-wave superconductor, and x<&$, where < x < i, for a d + is-wave 
superconductor. It should be noted that the minimum flux of x$q in the ground state for the d + is state occurs 
for the flux in a particular direction; the minimum flux in the other direction is (1 — x)<i>o, this asymmetry being a 
reflection of the broken time-reversal symmetry. (It has been noted previously in Ref. g| that the d + is state would 
lead to a ground state flux in a ring of x$o where x is neither nor i.) If the parameter L{Ja + 7b)27t/ (3>oc) is not 
large, the ground state flux through the ring will be less than i<f>o for a d X 2_ y 2 superconductor, as in the case of a 
single-junction ringoS, so that care will have to be taken to distinguish experimentally between this situation, and 
the one where the departure of the ground state flux from is due to broken time-reversal symmetry. These results, 
for the geometry of Fig. |H[ depend on symmetry arguments only, and do not rely for their validity on elementary 
models of the angular dependence of the Josephson critical currents. 

Finally, note that although the observation of a ring containing flux i n the proposed geometry would definitively 
rule out s-wave symmetry for the superconductivity, and would be consistent with d x - y 2-w&ve superconductivity, it 
would not definitively establish the superconductivity as being d x 2_ y 2-w&ve; d xy superconductivity would also be 



consistent with such a result, 
establish a eL.2_,,2 state. 



Other arguments, analogous to those of the Appendix, could however definitively 



VIII. CONCLUSIONS 



Since YBCO is orthorhombic, the order parameter describing its superconductivity is expected to transform 
acc<ffidi^[-Jj^WMj-af^t| ti ti irr educible representations of the orthorhombic point group. The experimental observa- 

nonzero bilinear Josephson critical currents along at least two of the three principal 
orthorhombic axes in tunneling from YBCO to a conventional s-wave superconductor is sufficient to identify the 
symmetry of the order parameter as being of the ux 2 + vy 2 (Ai g ) type (see Appendix). This provides the experi- 
mental justification for considering YBCO superconductivity as coming from a mixture of d x 2_ y 2-w&ve and s-wave 
contributions—which are the two tetragonal symmetry types compatible with orthorhombic A\ g symmetry. 

In a modeEj of mixed s and d x 2_ y 2 superconductivity in YBCO, there is a transition to a time-reversal invariant 
d + s state at the normal-to-superconducting transition temperature T c , followed by a possible transition to a time- 
reversal symmetry-breaking d + e l ^s state at a lower temperature. Corresponding to a d + e^s state there is the 
time-reversed d + e~ l( ^ s state of the same free energy density so that time-reversal domains may form in a sample. 
Since averaging over time-reversal domains (in the case where the size of these domains is much smaller than the 
size of the Josephson junction) makes the effects of broken time-reversal symmetry unobservable, the case of a single 
time-reversal domain is considered here. 
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Much of this paper studies Josephson tunneling from the d + e^s state to a conventional s-wave superconductor, 
and compares it with tunneling between a d + s state and an s-wave superconductor. The principal effect of broken 
time-reversal symmetry is that the Josephson junction acquires a characteristic phase which is neither zero nor tt. 
The differences in the characteristic phases of two junctions can be measured in a SQUID experiment. 

The symmetry characteristics of the Josephson currents in a twinned crystal can be understood in terms of the 
macroscopic invariance group of the superconducting state of the twinned crystal. For the case where the twins are in 



their time-reversal invariant d + s state (studied in Ref. 18) the macroscopic invariance group is that of a tetragonal 
d x 2_ y 2 superconductor for odd-symmetry twin boundaries, and it is that of a tetragonal s-wave superconductor for 
even-symmetry twin boundaries. Sincejua^Q^apic-d-wave behavior is observed in basal-plane Josephson tunneling 
experiments on-Jaighly twinned crystalsD'Qa'BBQ'BLlllj, the twin boundaries for the d + s state have to be of the odd- 
symmetry typeO, which implies that the c£-wave component of the order parameter is approximately continuous across 
twin boundaries. 

In the case of time-reversal symmetry-breaking d + e^s states, two types of twin boundaries are consistent with 
a continuous d-wave component, namely the odd-symmetry and the Ter twin boundaries. For odd-symmetry twin 
boundaries the macroscopic invariance group is that of a tetragonal d x i_ y i superconductor, but without the time- 
reversal operation. Detailed calculations confirm, that the Josephson tunneling for the d + s state and the d + e**s 
states are indistinguishable, and that they have the predicted macroscopic d-wave characteristics when the twins are 
present with equal weights. For Ta twin boundaries, on the other hand, the macroscopic invariance group is that of 
a tetragonal d + is superconductor. As expected from the macroscopic symmetry, in this case the c-axis Josephson 
currents cancel only partially when averaged over twins, and the Josephson junction characteristic phase difference 
measured in a corner SQUID experiment can differ from either zero or n. 

In order to determine the symmetry of the order parameter describing superconductivity in a single twin (or 
equivalently, the invariance group of this superconducting state), it is essential to do experiments on untwinned single 
crystals. Experiments on twinned crystals determine the macroscopic invariance group of the twinned crystal, which 
in general will be different from that of the superconducting state of an untwinned crystal, and which is determined 
by the symmetry of the twin boundaries. I— ■ 

A final section of the paper, stimulated by recent tricrystal-ring magnetometry experimentsc!3 on tetragonal 
Tl2Ba2Cu06+a, proposes and analyzes a ring geometry which allows s and d x 2_ y z superconductivity to be distin- 
guished by symmetry arguments only, and which makes unnecessary any assumption that an s-wave superconductor is 
unlikely to produce a characteristic phase of tt in a Josephson junction. The characteristics of d + is superconductivity 
for this geometry are also identified. 

ACKNOWLEDGMENTS 

This work was supported by the Natural Sciences and Engineering Council of Canada. Useful discussions with 
R. C. Dynes, A. E. Jacobs, and C. C. Tsuei are also acknowledged. 

APPENDIX A: ORDER PARAMETER SYMMETRY IN YBCO 

On the basis of experimental evidence, this section justifies the assumption that the superconducting cvrder. parame- 
ter in YBCO is a mixture of order parameters having d x 2_ y 2 and s wave symmetry. It is knownLJ'L3E3uj'E]EJ that the 
dc Josephson effect can be used to determine the symmetry of the superconducting order parameter. Since YBCO is 
orthorhombic, the logical first question is: Which of the irreducible representations of the orthorhombic point group 
does the order parameter belong to? The argument is straightforward and elementary but we have been unable to 
find it in detail elsewhere and therefore present it here. 

Consider a Josephson junction between YBCO and lead, the two superconductors being characterized by order 
parameters ipr and V'Pb, respectively. Here, ipr transforms according to the irreducible representation T of the 
orthorhombic point group £>2/i of YBCO and ?/;pb is assumed invariant under all transformations. Symmetry consid- 
erations show that the lowest order (in tp-p and V'Pb) contribution to the Josephson current across a junction will have 
the form 

Jn = iCl [M* Ph - # <Apb] , (Al) 
where n is the normal to the plane of the junction, and where J n is the current in direction n. 
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Table | shows which of the coefficients C£ are zero (by symmetry) and which can be nonzero, for each irreducible 
representation and for currents along each of the principal orthorhombic axes. Since only the ux 2 + v y 2 represen- 
tation in Table | can have nonzero Josephson currents along all three principal axes, and since nonzero Josephson 
currents haite-apparently been observed experimentally along all three principal axes in untwinned single crystal ex- 
perimentsueniJ, the ux 2 + vy 2 representation, also commonly called the A\ g representation, would appear to be the 
only acceptable one. 

When the bilinear contribution to the Josephson current is zero, however, the biquadratic contribution 



T (4)r • 



(A2) 



may be nonzero. Here, Cn is never required to be zero by symmetry. Close to the critical temperature T Cj pb — 7.2K 
for lead, the biquadratic critical current should have a temperature dependence j( 4 ) r cx |"0pb| 2 oc (T c .pb — T), 
i.e. linear in T c pb — T, whereas the bilinear critical current should have the square root temperature dependence 
J r oc IV'Pbl aiLm,- T) 1 / 2 . Measurements of the temperature dependence of the c-axis critical current by Dynes 
and coworkersEJOO clearly show the square root temperature dependence of the bilinear current. 

We are not aware of measurements which establish experimentally that the Josephson currents along the a and b 
axes are of the bilinear, as opposed to the biquadratic, type. An experiment which does this would be of interest. 

Note that the biquadratic Josephson current is proportional to sin[2(</>r — </>pb)], where (f>r — <f>pb is the phase 
difference between the two superconductors, whereas the bilinear Josephson current is proportional to sin(0p — </>pb), 
This difference means that a SQUID made from Josephson junctions for which the current is biquadratic would have 
a maximum current periodic in the flux with period <&o/2 rather than $o as is usually the case ($o is the elementary 
flux quantum). Also, the base frequency for the ac Josephson effect would be cjj = 4eV/fi rather than 2eV/h. For a 
microscopic calculation of the effect for tetragonal systems see Refs. H ||. 

If the bilinear character of the Josephson currents along the a and b axes of YBCO can be confirmed experimentally, 
then the ux 2 + vy 2 character of the superconducting order parameter in YBCO will be established. Of all the 
different irreducible representations of the tetragonal point group, only the d x 2_ y 2 and s representations are invariant 
under the operations of the orthorhombic point group. (A different way of seeing this is to note that ux 2 + vy 2 — 
d(x 2 — y 2 ) + s(x 2 + y 2 ), where the coefficients d and s can be found in terms of a and b.) Thus the ux 2 + vy 2 order 
parameter for YBCO can be considered to be a combination of d x 2_ y 2 and s (e.g. x 2 +y 2 ) basis vectors (and no other- 
type of tetragonal symmetry). These arguments assume a second order transition to the superconducting state. 



APPENDIX B: ANGULAR DEPENDENCE OF JOSEPHSON CURRENTS 

Figure ^| shows a Josephson junction between two identical tetragonal superconductors. The free energy per unit 
area of the Josephson junction in zero magnetic field has the form 

F = C(9,0')cos((f> -(/)'), (Bl) 

where the superconducting order parameters on the two sides of the junction are ip = \ip\e 1 ^ and ip' = \ip'\e % ^ , and 9 
and 6' define the orientations of the crystallographic shown in the figure. 

First assume that the superconductivity is either s-wave or d x 2_ y 2-w&ve. (Our definition of s-wave is that the order 
parameter is invariant with respect to all operations of the tetragonal point group; this is sometimes called gener- 
alized s-wave superconductivity. Similarly, our definition of d :E 2_j / 2-wave superconductivity encompasses generalized 
<i :E 2„j / 2-wave superconductivity.) Then symmetry considerations show that the coefficients must satisfy the following 
conditions: 

C(0, 6 1 ) = C(0', 9) = C{-9, -9') = C{9 + 7T, ff). (B2) 
Furthermore, for s-wave superconductivity 

C s (6 + n/2,6 , ) = C s (6,e'), (B3) 

whereas for d x 2_ y 2 superconductivity 

C d (9 + n/2,9') = -C d (9,9'). (B4) 
The use of these constraints yields the following Fourier series representations for the coefficients: 
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C s (6, 9 1 ) = V [C 4n ,4n' cos(4n<9) cos(4n'6>') 



+SinA n , sin(4n6») sin(4n 6>')] (B5) 
for s-wave superconductivity, and 

C d (6, 9') = {C±n+2An>+2 cos[(4n + 2)9] cos[(4n' + 2)9'} 

+S 4 „ +2 ,4n'+2 sin[(4n + 2)5] sin[(4n' + 2)0']} (B6) 

for d x 2_ y 2-w&ve superconductivity. The sums over n and n! are over all positive integers and zero; also, the coefficients 
have the property C^.j = Cjj and Sij — Sj t i, where i,j are 4n, An' or An + 2, 4n' + 2. 
The first few terms of these expressions are 

C s (9, 9') = Co,o + C 4 , [cos(4#) + cos(40')] + ■ • • , (B7) 

and 

C d {9, 9') = C 2 , 2 cos(2(9) cos(2<9') + 5 2 , 2 sin(2<9) sin(20') + . . . . (B8) 

Note that the expression proposed by Sigrist and RiceE] for the d x i_ y i case, namely Csr(#, 9') = C2.2 cos(2#) cos(2#'), 
has the property C*sr((9, 7r/4) = 0, which is not required by symmetry and therefore cannot be expected to be true in 



general. The simplified expression given in Eq. (B8), which contains two constants 6*2,2 and 5*2,2, does not have this 
difficulty. Also, if 5*2,2 = ~C*2,2i the sines and cosines in our two-term expression Eq. flB8| ) can be summed to give 
Cd = C*2,2 cos[2(# + 9')}, which is the expression proposed by Tsuei et alEtHi to model what they call the maximum 
disorder limit of a grain boundary Josephson junction. 
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TABLE I. Josephson Currents for Orthorhombic Symmetry. The value of the critical Josephson current is given for each 
irreducible representation V and for each of the orthorhombic principal axes, a, b, and c. 
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FIG. 1. Unit cells in neighboring twins under the assumption that the twin boundary (dashed line) is a plane of reflection 
symmetry of the crystal. 

FIG. 2. Schematic representation of d x 2_ y 2- and s-wave order-parameter components near odd and even reflection-symmetry 
twin-boundaries. 

FIG. 3. (a) The structure of the interface between a region containing (110) twin boundaries and a region containing (IlO) 
twin boundaries, as described in Ref. The rectangles indicate the relative orientations of the orthorhombic a and b axes 
in the different twins, and the plus and minus signs indicate the relative sign of the superconducting order parameter in the 
different twins. The c axis is normal to the plane of the figure. Part (b) of the figure is obtained from part (a) by a rotation of 
7r/2 about the c-axis. 

FIG. 4. Corner SQUID geometry for basal-plane tunneling between YBCO and lead. A and B indicate basal-plane Josephson 
junctions on two perpendicular faces of the YBCO crystal, while I a and Ib denote the corresponding currents. 

FIG. 5. Proposed geometry for ring experiment to test for s-wave versus d x 2 _ y2 -wa,ve superconductivity in tetragonal systems. 
SI and S2 are tetragonal superconductors with their c-axis normal to the plane of the figure and their basal-plane unit cells 
rotated by an angle n/4 with respect to each other. A and B are Josephson junctions; the angle a is chosen so that 
a tt/4. 



FIG. 6. Josephson junction between two identical tetragonal superconductors. 
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